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Problems of applying the classical kinetic theory to tj,e growth of small 
droplets from the vapor are examined. A solution for the droplet growth equa- 
tion is derived which is based on the assumption of a diffusive Held extending 
to the drop suvface. The method accounts for partial thermal and mass accommo- 
dation at the interface and the kinetic limit to the mass and heat fluxes, and it 
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utilizing directional fluxes yields boundary values in terms of known parameters 
and a new Laplace transform integral. Different initial assumptions coupled with 
this method yield various solutions derived by earlier workers. The solution is 
applicable to drop sizes both larger and smaller than the mean free path. 


KEY WORDS 


IS. DISTRIBUTION STATEMENT 


Laplace transform integral 
Cloud droplet growth 
Aerosol droplets 


19 SECURITY CLASSIF.(a(*AI*M,>«t' 


U nclassi fied — U nil mited 






Unclassined 


20. SECURITY CLASSIF. (ol thU ,««•) 

Unclassified 


mm 


21 . NO. or PAOCS 22. PAICE 

49 NTIS 


MSrC* Form I III (Rw l>te«mb«r 1 • ' : > 


For ule by Natiooil Tachnicil Information !>arvica, Sprinffiald, Virfinia 1 ) 1 5 1 








ACKNOWLEDGMENTS 


The authors express their gratitude to Gwenevere Jasper and 
Blair Herren for proofing the text and computations. John Hallett was 
supported, in part, by grant ATM-7920399 from the National Science 
Foundation, Division of Atmospheric Sciences. 


ii 



TABLE OF CONTENTS 


Page 


I. INTRODUCTION 1 

II. SOLUTION FOR THE GROWTH EQUATION 2 

A. Initial Assutnptiona 2 

B. Vapor Density Boundary Condition.....' 4 

C. Temperature Boundary Condition 9 

D . Solution of the Diffusion Fields 12 

III. NUMERICAL EVALUATION 21 

IV. COMPARISON WITH EARLIER STUDIES 29 

V. CONCLUSIONS 34 

APPENDIX EVALUATION OF THE INTEGRAL, A(a/A) 35 


REFERENCES 


40 



LIST OF ILLUSTRATIONS 


Figure Title Page 

1 A (a/X) as a function of v-^/X) computed fn^n the 

third-order approximation given in the Appendix, 
equation (A-6) 22 

2 Sample evaluation of B, given by equation (37), as a 

function of radius for venous values of the sticking 
coefficient, 3 23 

3 Sample evaluation of C, given by equation (38), as a 

function of radius for various values of the accommoda- 
tion coefficient , a 24 

4 Radial growth rate for a pure water drop as a function 

of radius, shown for various values of environmental 
temperature. 25 

5 Radial growth rate for a pure water drop as a function 

of radius, shown for various values of environmental 
pressure 26 

6 Radial growth rate for a pure water drop as a function 

of radius, shown for various values of ambient satura- 
tion ratio 27 

7 Radial growth rate for a pure water drop as a func- 

tion of radius, shown for various values of the sticking 
coefficient, 3 28 


8 Error in the radial growth rate caused by assuming the 
mass flux is proportional to the vapor density gradient 
rather than the mole fraction gradient. Curve A is the 
baseline, T = 20°C, P = 1000 mb, S = 1.01, 3 = 0.035. 

Curve B shows T = 0**C; C shows P = 500 mb; D shows 

S = 0.20; and E shows 3 = 0. 10 30 

9 Absolute value of the percentage difference between the 
Fukuta and Walter solution and the current solution. 

Curve A, 3 = 0.1, T = 20°C; Curve B, 8 = 0.035, 

T = 20®C; Curve C, 3 = 0.035, T = 0.0®C. The curves 
are dashed below 1.0 pm because the Fukuta and Walter 
solution was not intended to apply in this regioa 31 


iv 


list of symbols 


a 

A(a/X) 

B 

c 

C 

C(R) 

D 

F 

K 

L 

m 

in 


M- 

n 

P 

P 

Q 

r 

R 

R 

R 


w 


sat 


a 


w 

s 

S 

t 

T(R) 

Ta 


u. 


drop radius 

integral defined by equation (8) 

defined by equation (28) 

mean molecvdar speed 

defined by equation (30) 

molar concentration of the moist gas 

specific heat of air 

diffusivity of water vapor in air 
heat flux 

thermal conductivity of air 
latent heat of vaporization 
mass of the drop 
mass of air molecule 

molecular weight of water 

integer 

environmental pressure 
saturation vapor pressure of water 

heat energy 

radial coordinate from origin on drop surface 
radial coordinate from origin at drop center 
gas constant for air 

gas constant for water vapor 

ratio of drop radius to mean free path 
saturation ratio 
time, integration dummy 
temperature (variable with position) 
temperature of the air at the drop surface 

temperature of the drop 

temperature of the environment far from drop 
internal energy of air molecule 


V 



UST OF SYMBOLS (Contlnutd) 


V most probable speed 

X mole fraction of water vapor in the air 

mole fraction of vapor at the drop surface 

x^ mole fraction of vapor far from the drop 

a accommodation coefficient (fraction of heat energy exchanged) 

$ sticking coefficient (fraction of vapor molecules that enter drop) 

A "jump length" (denotes position where kinetic free stream 

region meets diffusion controlled region) 

c numerical constant of order unity 

n viscosity of air 

6 polar angular coordinate, origin at drop surface 

X mean free path 

71 3.14159 

p(R) vapor density (variable with position) 

p vapor density at the drop surface 

p vapor density of moist air 

D 

density of liquid water 

p^ vapor density of environment far from drop 

PsatCT) saturation vapor density for the temperature implied 
4> azimuthal angular coordinate, origin at drop surface 

't' mass flux of water vapor 

i^i(s,n) defined by equation (A-2) in the Appendix 


Vi 



TECHNICAL MEMORANDUM 


AN EX1ENDED CLASSICAL SOLUTION OF IHE 
DROPLET GROWin PROBLEM 

I. INTRODUCTION 


The growth or evaporation of water drops under motionless, steady 
state conditions in air has been the subject of a number of theoretical 
studies over the last 100 years. The problem is clearly of importance in 
cloud physics for the growth of cloud drops up to 10 ym and for the 
development or dissipation of fog and haze in the atmosphere. Similar 
applications can be made to aerosol/gas phase reactions and to numerous 
chemical engineering processes, the dispersion of agricultural sprays, 
spray drying, combustion, humidification, and spray cooling. 

The rate of growth or evaporation is limited by transport processes — 
of heat by conduction or vapor by diffusion in any inert gas present. 

The first solution to this problem is attributed to Maxwell [1]. He 
assumed that the surface of the drop produced a vapor pressure related 
directly to its temperature through the Clausius Clapeyron equation and 
that heat and moisture were transported by conduction and diffusion, 
respectively, in a spherical diffusion field from the drop surface. This 
first approximation is commonly quoted for growth of particles in cloud 
physics texts (2-5] because it provides sufficient approximation for par- 
ticles larger than approximately 10 ym. It was first noted by Langmuir 
[6] that these assumptions are not accurate, especially for small 
sizes. Three physical factors impact the boundary conditions which were 
used to represent the drop surface. First, there is not complete 
"accommodation” at the drop surface — not all vapor molecules which strike 
the surface stick to the drop, and the energy exchange when an air mole- 
cule strikes the surface may be incomplete. Second, there is a kinetic 
limit to the diffusional heat and mass transfer processes. These fluxes 
will increase as the gradients increase only to a point, the maximum kinetic 
free stream limit. Third, the drop disturbs the gas near its surface so 
that the "equilibrium” and "isotropic" conditions no longer strictly apply — 
the diffusivity and thernial conductivity are no longer constant and the 
velocity distributions are no longer Maxwellian. Considerable attention 
has been given to these difficulties. Many investigators have recognized 
that an exact solution is, in principle, not to be obtained using kinetic 
theory based on Maxwellian distributions and have sought a basic solution 
using Chapman Enskog theory, neutron scattering theory, or similar 
approaches to the fundamental kinetics. This work has yielded valuable 
insights into the problem, but it has not yet produced a complete solution 
to this difficult problem. With the exception of the first-order attempt of 
Monchick and Reiss [7], this work is beyond the scope of this study. 
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Other investigators [8-18] have developed solutions using classical kinetic 
theory methods coupled with various eqpproximations and semiempirlcal rela- 
tions to yield solutions with slight differences in empiricism and level of 
approximation. 

The primary problem here is that there is seldom a clear picture 
of why these differences exist, since the classical physical assumptions 
which form the basis of these approaches are very similar or even identi- 
cal from author to author. Thus the question remains as to which solution 
most accurately reflects the physical situation. Much of the problem 
results from the usw of an artificial computational device, the "jump 
length," A , to reach a solution in analogy with the problem of a hot fila- 
ment in a rarefied gas. This approach divides the environment around 
the drop into two regions, a near zone where the vapor flux is controlled 
by free stream kinetics and a far zone where the flux is controlled by 
diffusion. The "jiunp length" is the di'^tance from the drop surface to 
the boundary between the zones. The jump length must be evaluated to 
reach a solution of the drop growth problem, but no clear line of physical 
reasoning has been presented that leads to a definite relationship to the 
known parameters: the length appears to be only a computational artifice. 
Thus, various assumptions have been made, and a proliferation of solutions 
is the result. 

This report presents a new solution to the droplet growth problem 
that extends the solution to drop sizes smaller than the mean free path 
and that does not rely on the "jump length" artifice. The solution is 
based on the same classical physics and a boundary value method which 
closely parallels the original Maxwell method. The impact of the various 
assumptions can be clearly seen in the Hnal solution. Thus, the method 
is useful for pedagogical reasons. More importantly, however, it is shown 
how other initial assumptions yield several of the more popular solutions 
derived using the "jump length" approach. Thus, the physical basis for 
some of the discrepancies currently in the literature is made clear. 


11. SOLUTION FOR THE GROWTH EQUATION 


A . Initial Assumptions 

Consider the growth of a drop of pure water (radius = a) at rest 
in an environment of air and water vapor. The well-known correction for 
a dissolved impurity in the drop could be easily included, but it has been 
omitted here for the sake of clarity and brevity. Note that since a clas- 
sical kinetic theory solution to the problem is sought, the accuracy of 
this treatment, like any heat or mass transport problem, is limited to a 
few percent at best because the temperature and density gradients which 
produce the transport also invalidate the basic assumptions upon which 
the theory rests. Thus, as a general guideline for this study, approxi- 
mations and assumptions which could alter the solution by less than 
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1 percent have been freely made, while larger factors have been treated 
as rigorously as possible. The resulting composite theory is expected to 
be accurate to within approximately 3 percent for all drop sizes. 

Consider now the three physical factors mentioned previously which 
affect the Maxwell-type solution. The first two factors. Incomplete accom- 
modation and the kinetic limit, represent cutoffs which keep the flux den- 
sities from becoming infinite at small drop sizes. Fractional sticking of 
molecules and exchange of energy proportionally require greater fluxes 
to obtain the same result. Thus, the accommodation coefficients enter as 
multipliers in the theory and their consideration is clearly important unless 
they are known to be near unity. The kinetic limit becomes important 
for drop sizes approaching the length of the mean free path in the gas, 

X, as can be seen readily by evaluating when the kinetic free stream flux 
equals the diffusion flux. The diffusion flux is approximately proportional 
to the vapor density gradient 


cAp r»Ap _cXAp 

T' " ^ ^ “ T IR • 


Here we have made use of a classical kinetic theory^ relationship between 
the diffusivity, D, and the mean molecular speed, c [19]. Thus, 

AR - 4X/3. Since we seek a solution applicable for small drop sizes, this 
factor must also be included in the analysis. 


The third factor to be considered is the change in the gas proper- 
ties near the drop which result from the composition and temperature gra- 
dients in that region. There are first -ojrder effects of these gradients — 
heat conduction and vapor diffusion — which obviously must be included 
in the theory. However, these gradients also cause the thermal conduc - 
tivity and the vapor diffusivity to be functions of distance from the drop, 
they invalidate the Maxwellian molecular velocity distribution in that 
region, and they give rise to "thermal diffusion" of the vapor (Soret 
effect) and a heat flux due to the concentrati n gradient (Dufour effect). 
All of these are second-order effects which make a contribution that is 
small compared to the errors caused by the basic limitations of the clas- 
sical kinetic theory. For example, examination of an extreme case for a 
drop in the atmosphere — growth of a pure water drop by the Maxwell 
equation under a constant 10 percent sapersaturation starting from its 
critical radius, 0.01 ym — shows that the temperature deviates from the 
environmental value by only a few tenths of a percent and the water 
vapor concentration, a minor constituent in the gas, deviates by only 


3 percent. Because the diffusivity varies as the where is 

the gas density , and the conductivity varies approximately as /Y, 
assuming that these quantities are constants, will not contribute appreci- 
ably to the error already present in the classical approach. This allows 
us to greatly simplify the problem because it implies that Laplace’s equa- 
tion for the temperature and vapor density is satisfied throughout the 
region surrounding the drop. 
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In summary, the following assunsptlonB are made in this analysis: 

1) Quasistatic equlUbiium: i.e. , changes in drop raditut and other 
conditions occur so slowly that equilibrium kinetics can be applied. For 
justification of this assumption see Reference 2. 

2) The air and water vapor behave as inert, ideal gases. 

3) Maxwellian kinetic theory applies. 

4) The thermal conductivity and vapor diffusivity are assumed to 
be constants, independent of distance from the drop. 

5) Pick's law of diffusion and the heat conduction equation apply 
throughout the environment of the drop, even near the surface. Thus, 
the only discontinuity is at the drop -air interface. 

6) The drop temperature T^, the mean temperature near the drop 
T, and the environmental temperature T^ are sufficiently close that 


where this factor appears as a multiplier. 


B. Vapor Density Boundary Condition 

Far from the water drop the water vapor density has a constant, 
known value, p^. As one approaches the drop along a radial vector, the 

vapor density must be a smooth, continuous function of radius, R, until 
the water surface is reached since there are no sources or sinks outside 
the drop and the environment is assumed unifonn (constant diffu ivity). 
Approached in this way, the limiting value of p(R) at the drop surface 
is defined as p . p does not equal the equilibrium value of vapor density 

over the drop surface, as is assumed in the Maxwell-type solution. 

The difference is the "jump" which is well known from kinetic theory. 

Note that this is a jump in vapor density at the location of a physical dis- 
continuity in the mMium, the drop surface. It is not e jump in position, 
p is the boundary condition for the diffusion field. It must be solved 

A 

for in terms of known parameters and parameters that can be eliminated 
in the subsequent derivation. 

Consider first a drop in equilibrium, dm/dt ^ c. From elementary 
kinetic theory one can write the water vapor mass flux (per unit area 
per unit time) against the drop. 
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( 1 ) 
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(sat) _ 
in 


^aat 




because in this case the vapor density is a constant, Equation (1) 

holds subject only to the static equilibrium assumption, independent of 
the sise of the drop. 

An expression for the mass balance at the air-droo interface is 
also easily obtained for the equilibrium drop: 


Flux in s Flux out = Flux reflected Flux leaving drop 


Making the definition, 3 = the fraction of incident water molecules absorbed 
into the drop, it follows that 


= a - 


3 ) 


(sat) 


in 


+ t 


Thus, for the flux leaving the drop 


= 3 


(sat) 


in 


( 2 ) 


For a drop not in equilibrium, i.e. , growing or evaporating, the 
mass balance condition at the interface is given by a similar equation. 


Flux in = Flux reflected + Flux leaving drop 
+ Rate of mass change /area; 


in 


= (1 - 3 ) 4 


in 


+ >}> 


4t«' 


dm 

dt 


Here it can be assumed that the flux leaving the drop, <1’^, is 

only a function of drop temperature and not a direct function of the equi- 
librium condition of the drop. Thus, one may employ equation (2) and 
arrange terms to yield 
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1 


(3) 


♦ ^ ^ <*«'> + 
*in 


4nv 


dm 

j ar 


t 


where is to be computed with equation (1) using the drop 

temperature. 


Next an equation is required that is similar in form to equation (1) 
to relate the flux in to the vapor density field. Equation (1) no longer 
e^plies because p(R) is not constant near the drop. Thus, it is neces-- 
sary to return to one of the integral formulations from which equation (1) 
is normally derived (Loeb [ 19J ) and make t*^ a adjustment for a vapor den- 
sity which is a function of R 


, - . ,3 -c*/,2 ■ 


! 

c=0 


d. { j j 

r =0 »=0 < 1=0 


^1* t A 

o(R) e 8ln0 oo?0 d<>d0dr 


(4) 


=» n/2 2fT 

^in = TiX / / / 

r-O &=0 4>=0 


p(R) e sin0 cos0 d(f d0 dr . 


The integration is over a half space in a (r, 6, 0) coordinate sys- 
tem centered on the surface of the drop. Here X is the mean free path 
in the gas, v is the most probal,le speed, and c is the mean molecular 
speed. From elementary kinetic theory, 


c = /sr^T/tt 


(5) 
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In Section II-D it is shown that Pick's Law for the diffusion field 
implies that the vapor density obeys the relation 


P(R) = P„ 


+ <Pa ■ P«) «/R + [Pg “ P(R)1 



a/R . 


Sines I - p(R) I < I Pg - P„ I R ®® noted nreviously, 

fT„ - T )/T is less than a few tenths of a percent in even extreme 
cases,* the last term can be ignored; thus 

p(R) = P„ + <Pft " P«,^ 

Using the law of cosines, 


R = (a^ + + 2ar cos0)^^^ 

Performing the 41 integral, breaking the integrand, and performing the 
elementary c and 0 integrals in the first term yields 


Pco 3 ^Pa "Poo) c A (a/A) 

‘t’ . = — }j — + ft 

in 4 2 


We have made the definition 


(7) 


it/2 

A(a/A) s a/A / / 

0 0 


e ' sin9 cos0 d0dr 
(a^ + r^ + 2ar cos9)^^^ 


( 8 ) 



• For droplet growth. The error could reach 10 percent for drops 
evaporating in very dry air. 



where E 


i(*)’ { 


-at/x 


dt 


The expression (9) represents the second-order approximation to 
the value of the integral. It 0ves a value for A(a/X) which is 10 percent 
too low for a s X, but the accuracy improves rapidly for larger values of 
drop radius. Note that an error in the value of A(a/X) causes a much 
smaller error in the drop growth rate. Evaluation of the integral and 
expressions for successively better approximate values and upper and 
lower bounds are discussed in the Appendix. 

Substituting equations (7) and (1) into equation (3), one approaches 
the final solution for the boundary condition. 


c (p - p^) cA(a/X) 

4 - ■■ . 

4 2 



( 10 ) 

Using equation (5) and solving for p , making the approximation 


gives. 



1 \ 1 / 2ti dm 

2A(“)j S.Ka^a) / dt ' 


( 11 ) 


Equation (11) represents the boundary condition for the vapor 
diffusion field at the water-air interface. Notice that in the limit of 
large drop radius A(a/\) + 1/2, so the boundary condition becomes 





va » X) 


( 12 ) 
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C. Temperature Boundary Condition 

The solution for the temperature boundary conditions very closely 
parallels the vapor density case. If a is defined as the fraction of energy 
exchange for the air striking the drop surface » 


- Actual energy exchange 
“ “ Mai^mum pos^le exchange 


consideration of the energy flux balance in equilibrium yields, in analogy 
to equation (2), 


^d = 


(13) 


for the flux leaving the di p. Similarly, for a growing or evaporating 
drop the flux balance condition is 


^in 


= (1 - a) 


*^in 


+ - 


5 

47ia^ 


dm 

dt 


where L is the latent heat of condensation. Assuming the directional 
energy flux emerging from the drop is independent of the state of equi- 
librium, one obtains 


^in = ^'in<««^> 


•iriaa 


dm 

dt 


(14) 


Now expressions must be obtained for the incident energy fluxes 
for the growing and equilibrium drops. In both cases we begin by noting 
that the energy flux, F, equals the product of the energy per molecule 
times the molecular flux. Using equation (4) for the molecular flux with 
the vapor density replaced by the gas density and writing the energy 
per molecule as translational energy plus iixternal energy to separate out 
the dependence on molecular speed, yields 


9 


, - n/2 I* - /m.c* iU, \ , .„*/„> P|,<") ../» 

^ 1 -Sr* 


sinO co«0 dcded^dr 


Integrating over c 


P - 1 p»<R) 


sine cose d(^d6dr 


Since c = 2v//F , = 2R T, and c„ = 3R„/2 + 3U./3T, where 

a V a a 1 

R is the gas constant for air, 

ol 


® it/2 2it 


F = - 


( V t £t ti n 


0 0 0 


T(R) P (R) e sinG cos0 d4>d0dr 

o 

(15) 


In equilibrium T(R) and p (R), the gas density, are not functions of 
position, so the elementary integration gives 


Fi„(sat) = (C^ , R^/2) Pg /rj^T^) 


' <3R.T^) /RJ^/(2.) , 


(16) 


where is the drop temperature (which equals the gas temperature). 
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For drops not in equilibrium the analysis of the diffusion field 
given in Section Il-D shows that the temperature field is given by 


T(R) - + (T^ - T„) a/R 


( 17 ) 


where is the gas tmnperature at the drop surface and is the envi- 

a OP 

ronment temperature. After inserting equation (17) into equation (IS), 
integration jdelds the desired solution for the energy flux into the drop. 
Fin> In this case the gas density is constant in the region about the 

drop (to a close approximation), but the temperature shows the same 
radial dependence as the water vapor density did in the previous case. 
The solution is 




SR c p 


81 


2n tt/2 r 1 r/\ 

J / 1 T„ + (T - T^) S ® sinecosfl drdUd^ 

0 0 0 a -x- R 

W m 


f3R,Pgl 


T c 

CO 


(T^ - T^) c A(a/\) 


(18) 


A(a/\> is again defined by equation (8). Combining equations (14), (16), 
and (18) as before, one obtains the final boundary condition: 




L 

24«c.Ra 



Again, since A(a/\) 


1/2 for a » A, the large drop limit is 


T 


a 


‘d 


L T / . V 

a / dm \ 

12Tiaa^R p \ ^ / 

« « 


(a » X) 


(20) 
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D. Solution of the Diffusion Fields 

Now that the boiindary conditions for the vapor and temperature 
fields have been established, equations (11) and (19), the solution for 
the droplet growth rate can be obtained. The problem contains four 
unknown parameters: the growth rate, the latent heat released at the 

drop, th'> drop temperature, and the equilibrium vapor density above 
the drop. Thus, four equations are required: Pick's law for the diffusion 

of vapor, Fourier's law for the conduction of heat, the proportionality 
between latent heat release and change in mass, and the Clausius- 
Clapeyron relationship between temperature and equilibrium vapor density. 
These yield a solution in terms of the known environmental parameters. 

The method parallels the textbook solution of Mason [31 with different 
boundary conditions and with Pick's law written in terms of the mole 
fraction, x, rather than the approximate expression in terms of vapor 
density. 

Beginning with Fourier's law of heat conduction and noting that , 
for a steady state system, the heat flux across any spherical shell con- 
centric with the drop is constant, independent of time and radius, we 
have 


^ = 4ttrV = -4itR^K ^ = constant. (21) 

dt dK 

An additional term could be included in equation (21) to account 
for the energy flux carried by the mass flux. However, this term must 
be less than dm/dt, where is the specific heat. Since the total 

heat flux must equal L dm/dt, the fraction contributed by the mass flux 
term must be less than C^/L ~ 0.0005, which is clearly negligible. There- 
fore, returning to equation (21), 


19 f dR 

J r 2 


-4ttK 



Integration and application of the boundary conditions, T = 
<», and T = T at R = a, quickly yields 

cl 


at 

00 


^ = 4,»K (T^ - T„) . 


( 22 ) 
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and 


T = T. ♦ (T, - T.) I . 


which is equation (17) which was used earlier. If one makes use of the 
definition of latent heat, L, equation (22) can be related to the growth 
rate 


dm 4itaK 
dt ■ 


TT - T„) 


(23) 


Next a similar solution for the diffusion field is obtained using Pick's 
law (201 and the "steady state" assumption. 


4irR“<f = constant 


4,r2c(R)M„D 
(1 - X) dn 


(24) 


where x is the mole fraction of water vapor, C(R) is the molar concentra- 
tion of the moist gas (air + water vapor), and M is the molecular weight 

w 

of water. From the gas law , 


C(R) = P/(R^T(R)) . 

Inserting equation (13) for T(R) into (24) and integrating from 
R = a to «> gives 
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Or. .'^er integration. 





Solving for dm/dt and neglecting the higher order terms leaves 


dm _ 8iraPD /*■» ~ *e \ 

ar = (T. + T,) ^ I - X. ; 


(25) 


This approximation is quite good. The series expansion is always valid 
since the absolute value of (x_ - x„)/(l - x ) is less than one. In the 

rather extreme case of a drop at 20®C evaporating into dry air; 
x = 0, x = 0.023; the error caused by neglecting the second-order term 

is only 1.2 percent. 

The mole fraction is related to the vapor density by the expressions. 


’^a = 


^w^a^a 


= 




Therefore, 


. SiraD (P^T^ - p„T„) 

am 00 (» a 

dt + Tj,) (1 - R^P^T„/P) 


471sD 


( t . [ 


<0. - o.) - r 

OO 


] 


41^0 


<0. - f=a> - r ‘■'a ■ 

OO 


(26) 
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Ibe aiJproKlmtloQs used to reach the fiaal fom of equation (26) 
introduoe negligible «*ror ae long as the quaeletatlo assuniition is 
valid. Ihey can be dioen to be equivalent to neglecting the tern 


(P^ - P.) (T^ - T.)/ 2T. 


plus higher order terns. Considering a dn^ at 20‘*C evaporating in 
90 percent relative hunidity - a significant departure from quasi- 
static - the ezTor oontributed by this approxiinatiGn would be less than 
0.7 percMt. 


Because of the spherical symmetry and steady stati assumptions, 
the mass flux across any spherical shell concentric with the drop is a 
constant. Thus, rewriting equation (26) for i^’hitrary R leads to an 
expression for the vapor density profile. 


dm 

= constant 


” [«>- ‘ - 17 ■ 'f-'] 

^ >T<H) - T.l] . 

Solving for p(R), using equation (17), gives 
P(R) = P« + <Pa ■ 

~ p + (p - p ) 


1+ [pg - P(R)1 ( ”T~/ R 

a 

R 


IS 



Thia is the vapor density profile, equation (6), which was employed 
earlier in Section II-B. 

Now p can be ettminated from equation (26) by substituting in the 
boundary condition, equation (11), 


dm 4 toD 

3T - / 




Next, substituting from equation (23), (T - T ) is eliminated to give a 
solution for dm/dt. ® “ 


dm 

HT 


4™D (p„ - 




^ur Pc»*^ao 

1- 


D^r/R^T„ 


28^ 




DP„L 

*icr- 


4-naD 

B 


- Paat> 


(27) 


where the parameter B is defined to equal the denominator of (27) 



D,^^^r/R DLp„ 

Yi oo 00 Yf CO 00 

p iBaA(a/\) * TTT7 


(28) 


This expression, equation (27), appears in the original solution by 
Maxwell [1] with B ^ 1. It represents a solution for the droplet growth 
rate which is independent of the thermal conduction field. Thus, one 
must have independent knowledge of the drop temperature to specify 
Pgat equation (27). 

Likewise, an expression for the growth rate can be obtained from 
the thermal field by substituting the temperature boundary condition, 
equation (19), into equation (23). 
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( 29 ) 



27TKa (T^ - T^) 



* 6a«R,p^f) J 


" LC 


T«) 


» 


where, to simplify notation, we have defined 



Just as the saturation vapor density over the drop, Pg^^t must be 

known to use equation (27), the drop's temperature must be known to 
use equation (29). This difficulty can be eliroinatsd by using an expres- 
sion which will be derived from the Clausius-Clapeyron equation. 



dP 


sat 


dT 


R T' 
w 


(31) 


Using the gas law to wriie this in terms of vapor density, rather 
than pressure, gives 


^Osat'T) 



Integrating from to yields: 


In 


'Oa,t"d>\ '-‘I'd-T.) 



(32) 
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The approximation and the expansion, ln(x) - x - 1, 

for X s 1 applied to equation (32) yields 


(I- - twT.) (Tj - T. 


Taking the exponential of both sides, expanding, and dropping the 
higher order terms leaves; 



(T^) 

(T„) 



Substituting for (T^ ~ equation (29) leaves 


‘^sat^'^d^ ~ 


^sat^^< 




~ ^w^oo \ / LC \ /dm \ 

■ \h„t! 


(33) 


Recalling equation (27) and dividing it by Psa^(T^). one has 


^sat _ B /^•"\ 

Vdty 


(34) 


Since it is assumed here that the water drops are pure, the equi- 
librium vapor density over the drop, Pg^^. is related to the equilibrium 

vapor density over a flat surface of water at T . , p ,(Tj), by the Kelvin 

Q Sflf Q 

equation: 


'’sat = *’“’(.'LR°T^a) 
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* p, R T a 
L w <» 


( 35 ) 


The linear approximation for the expon^tial is adequate for wacer 
drops larger than O.I ym. For aj^lication to smaller drops the ourjdratic 
teim should also be included. The teim for dissolved ui|)urlties could 
have been Included at this point by the textbook method [3,4], but it is 
omitted here to ^K>rten the solution. 

Substituting equation (35) into equation (34) and adding result 
to equation (33) gives: 


*^sat' p, R T a 

L w °° 

Psat^*"-^ 


* (4«aDp^t<’'„)) (ar ) 


Using the definition of the saturation ratio, S, 


Psat^'*^oc> 


The final solution can now be written: 


dm 

4ira 

(S - 1) - 

2o 


dt ■ " 

’lc/' 

- 

L / B 



R t2 j 




W CX> ' 

\ i 


r 


(36) 
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where 


B 


. u(^) [l - 


d/2it/R^T„ DLp„ 

^ ~m ATa AT ^ ■RT” 


P T 

'^OO 00 


( 37 ) 


/axf 

C = 2A^^) 1 + 


/2ir/R^ 


gA(aA) J 


( 38 ) 


A(a/X) is the integral, equation (8), which is evaluated in the Appendix 


/•“ sine COS0 dr de 

A(a/X) - a/A / / - -r ^ 

0 0 (a + r^ + 2ar cose)^'^ 


Since 


4 3 

in — ^ TT a p , 


dm/dt = 4ita Pj^ da/dt , 


the radial growth rate can be written 


a 


dt 


(S - D- 


20 


PlLC /L- yrTY 
^ I R T ^ J 

\ w “■ / 


PlB 


( 39 ) 


As noted previously, a second-order correction to the Kelvin tern may be 
desirable for sane i^lications involving very smell droplets. 
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III. NUMERICAL EVALUATION 


The numerical values of A(a/X), B, C. and da/dt are illustrated 
in Figures 1 through 7 to indicate their magnitudes and dependence upon 
temperature, pressure, saturation ratio and drop size. Figure 1 illus- 
trates the value of the integral, A(a/X), as a function of drop radius 
divided by the mean free path. The third -order approximation given in 
the Appendix, equation (A-6), was used to compute the graph. Note 
that the value of the integral is well defined for all values of radius. 

In the limits, A(a/\) 0 as a 0 and to 1/2 as a | 

The difference between this formulation of the droplet growth theory 
and earlier solutions is contained in the terms "B" and "C” which are 
given by equa’ >ns (37) and (38). These are plotted against radius in 
Figures 2 and 3, respectively, for various values of the thermal accom- 
modation coefficient, a, and sticking coefficient, g. If a and g were 
about equal, then "B” and "C" would be of nearly equal significance; but 
since it is generally believed that a is near unity and 6 is small , probably 
between 0.03 and 0.04, the "B” term must be regarded as having the 
greatest impact on the growth rate. However, as pointed out by Chodes 
et al. [21] and others, the measurements upon which this judgement is 
based are difficult to make with high accuracy , so the values of a and 
3 must still be regarded ;viih some suspicion. In the limits both "B" 
and "C" become inversely proportional to "a” for small drop sizes; but, 
since the entire denominator of the growth rate equation is multiplied by 
the radius, the growth rate would remain bounded as a -»■ 0 except for 
the Kelvin term in the numerator. As the drops become large, ”B" 
(approximately) and "o" approach unity; thus, equation (39) approaches 
the c.ussical solution. 

The remaining figures (4 through 7) illustrate the magnitude of the 
radial growth rate, equation (39), as a function of radius with tenperature, 
pressure, saturation ratio and sticking coefficient as parameters. As 
will be discussed in the following section, the numerical values do not 
differ appreciably from those presented by other authors. In each of 
the curves plotted, the growth rate is carried to zero by the Kelvin 
term. If that were absent, the curves would approach a constant value 
somewhat greater than the peak values shown. Of course, to accurately 
reflect the growth of a real drop in the atmosphere one must include, in 
addition to the Kelvin term, a term to describe the influence of the nucleus 
and any chemical impurities in the drop. 
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(oA) 

Figure 1. A(a/A) as a function of (a/X) computed from the third-order 
approximation given in the Appendix, equation (A-6). 
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RADIUS 0<m) 

Figu!*e 4. Radial growth rate for a pure water drop as a function of radius, shown 
for various values of environmental temperature. 




RADIUS 0<m) 

Figure 5. Radial growth rate for -i pure water drop as a function of radius 
shown for various values of environmental pressure. 




VARIOUS SATURATION RATIOS 
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RADIUS (Mm) 

Figure 6. Radial growth rate for a pure water drop as a function of radius, shown 

for various values of ambient saturation ratio. 




FOR VARIOUS SnCKiNG COEF 
P«1 ATMOS. 



RADIUS (fiml 

Figure 7. Radial growth rate for a pure water drop as a function of radius, shown for 

various values of the sticking coefficient , 3 . 




IV. COMPARISON WITH EARLIER STUDIES 


Sorting out the differences in the notation, assumptions, and lines 
of reasoning which were used in the various prior treatments of this 
problem is a difficult and generally profitless enterprise because the dif- 
ferences are often quite subtle and not verifiable by experiment. Thus, 
a comparison of this type will not be attempted in this report. However, 
it was noticed in the course of this research that several of the more 
recent and widely used expressions for the droplet gp:«wth rate could 
easily be derived using the line of reasoning presented here — with 
appropriate changes in the initial assumptions. Thus, it becomes an easy 
matter to evaluate the differences among these formulations. It is these 
cases which are discussed in this section. 

Consider first the solution presented by Fukuta and Walter in 1970 
[14]. It differs from the solution presented here in two aspects: (1) the 
terms - R^^^PooTj^/P and DLp^/KT^ in our expression (37) for ”B” and 

(2) A(a/X) = 1/2 rather than the integral, equation (8), presented here. 
The additional two terms in the expression for "B” in the solution pre- 
sented here derive from the use of equation (24), the basic formulation 
of Pick's law of diffusion in terms of the mole fraction rather than the 
approximate formulation in terms of the vapor density gradient. The 
propriety of using the mole fraction formulation was pointed out by Peard 
and Pruppacher [22] and Carstens [18]; its influence on the growth rate 
is illustrated by Figure 8. The second condition, substituting 0.5 for the 
value of the integral is equivalent to using the large radius limit. As 
Figure 1 illustrates, A(a/X) approaches to within 10 percent of this value 
by the time a/X reaches six. Since the Fukuta and Walter solution is 
limited to a >> X by their method of derivation, the solution presented 
here can be viewed as an extension of their work to small drop sizes. 
However, it is perhaps more revealing to note that making the assumption 
that the vapor density and temperature are constant near the drop, i.e. , 
p(R) = p in place of equation (6) and T(R) = T in place of equation 

(13), yields A (a/X) = 1/2 directly. Thus, the Fukuta and Walter solu- 
tion fails to account for the vapor density gradient near the drop which, 
it turns out, is not important unless a < 6X. Figure 9 illustrates the 
absolute value of the total difference between the two solutions. 

Another solution which is closely related to the Fukuta and Walter 
formulation is the one developed by Rooth [ 12] . Rooth's solution is equiv- 
Eilent to equation (27) if one seis A(a/X) = 1/2 and drops the two terms 
discussed in the preceding paragraph. Rooth considered only mass trans- 
port and ignored the heat transport portion of the problem. Thus, his 
solution and others similar to it can be applied when the drop temperature 
is known independently, as when the drops are resting on a surface or 
suspended from a thermocouple. As far as the treatment of the mass 
transport problem is concerned, the Rooth solution is equivalent to the 
Fukuta and Walter treatment: therefore, the same comments apply to both. 
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RADIUS Omi) 

Pigrure 8. Error ir the radial growth rate caused by assuming the mass flux is proportional to the 
vapor density gradient rather than the mole fraction gradient. Curve A is the baseline, 

T = 20°C, P = 1000 mb, S = 1.01, 3 = 0.035. Curve B shows T = 0®C; 

C shows P = 500 mb: D shows S = 0.20; and E shows 3 = 0.10. 





A third formulation which is widely used and is also closely related 
to the Fukuta and Walter work was developed by Carstens [17,18]. 
Carstens employed a classical kinetic theory approach from Kennard [23] 
to account for the temperature and vapor density gradient near the drop. 
The method presented here yields the Carstens solution if one replaces 
equation (6), the solution of Laplace's equation near the drop, with an 
expression of the form 


p(R) = p + R CO80 


dp 

dR 


(40) 


and likewise for the temperature field. This is equivalent to considering 
a flat water surface and a one-dimensional variation in vapor density when 
deriving the boundary conditions. Thus, Carstens fails to correctly 
account for the spherical geometry near the drop. In his second paper 
[18] he works from the mole fraction formulation of Pick's law, but he 
makes approximations so the the "mole fraction" terms do not appear in 
his equivalent of our "B", equation (37'. The droplet growth solution 
developed by Monchick and Reiss [7] bears the same relationship to the 
Carstens solution as the Rooth solution does to Fukuta and Walter. Thus, 
the preceding comments apply to it as well. 

The difference between Carstens* work and the solution presented 
here can be evaluated by considering the second-order approximation to 
the value of A(a/A) presented in the Appendix, equation (A-5). For 
a/ A >3 the three exponential and exponential integfral terms provide a 
negligible contribution and one has 


A(a/A) ~ 1/2 - A/3a . (41) 


Making use of the classical kinetic expressions D = cA/3 and 
K = I enR„ where n is the viscosity (1/3 pjsX) and e is a nurneiical con- 

4 a g 

stant between 1 and 2.57, depending upon the mathematical analysis (see 
Loeb, [19] pp. 245, 441), one quickly obi ins 

Aj = (3D/4) /IT/R^ (42) 

and 


X - 3K 
2 " 10 


/ 27T/R T 
a 


(43) 
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Thus, one can write 


A(a/A) - 1/2 - 4D/a / 2ir/R^T 


( 44 ) 


and, choosing e ~ 1.2, 


‘ - nsfbr 



(.5) 


Substituted into equations (37) and (38), respectively, these equa- 
tions yield the Carstens solution. Equation (42) is an expression for the 
mean free path of the water vapor molecules, while equation (43) refers 
to the air. 

Relations (42) and (43) ca.i also be used to aid intercomparisons of 
the solution presented here with those developed by Fuchs [9], Bradley 
et al. [10] and others that leave the final expression written in terms of 
X or A , the "jump distance" which is not deHned analytically in the solu- 
tion. Consider, for example, the well-known Fuchs solution [9]. Fuchs 
considers the vapor and heat transport problems separately and does not 
combine the two, so we compare with our expression (27). Fuchs* expres- 
sion (5.6) in our notation is 


dm 

ar 


\a + A 6a / R^T J 


47TaD(Pg^j - p^) 


2^F ^ eeAp /^) 


if one defines Ap = 1/2 - ^ , which compares to our A(a/X). Fuchs 

argues, based on analogy with heated filament work, that A = tX where 
£ is a number greater than one for small drops. This does not appear 
correct in light of the present analysis which shows that Laplace's equa- 
tion for the thermal field is obeyed to within a few percent even for 
rapidly growing or evaporating drops. The problems associated with the 
sharp temperature differences near a hot filament do not enter the analy- 
sis. Fuchs notes that > 2/3 for large drops, so in this case 


33 



A„ ^ 1/2 - X/3a = 1/2 - 4D/a / 2 tt/R.„T 
r w 


This is, again, equation (44). Thus, for large drops the Fuchs 
solution reduces to the Carstens. Monchick and Reiss solution. As dis- 
cussed previously, it differs in the presence of the two terms related to 
the use of the mole fraction gradient form of Pick's Law and in the value 
of A(a/A). 


V. CONCLUSIONS 


In this report a solution for the g^wih of water droplets under 
quasistatic conditions is presented which is accurate to within the basic 
limitations of classical kinetic theory, approximately 3 percent. The anal- 
ysis is based on the concept that the diffusion field and thermal conduc- 
tion field are continuous right to the drop surface; there are no abrupt 
changes in the properties of the media in the near vicinity of the drop. 
The validity of this concept is verified by the analysis for essentially all 
circumstances of interest in the atmosphere; that is, that "quasistatic" 
approximation implies that the temperature and vapor density fields near 
the drop are Laplacian to within an accuracy of just over 1 percent or 
better. For a drop evaporating in a very diy environment, the vapor 
field is non-Lapladan to an amount approximately equal to the dew point 
depression divided by the absolute temperature, i.e. , a figure which could 
reach 10 percent. 

The theory presented here is rigorous within the limitations of the 
classical kinetic theory. Thus, no approximations have been made which 
would cause a deterioration in the accuracy beyond the few percent level. 
The solution is applicable to drops small compared to the mean free path 
as well as to laige drops; thus, it represents an improvement over the 
previous theories which required that the drop radii be large compared 
to A. 


Numerical evaluation and study of the differences which exist 
between this and previous theorie:. reveal that the differences are not 
significant in most cases. Indeed, since it appears that the state of 
the art limits the precision of g^rowth rate measurements of small diops 
to something in excess of 10 percent and that a and t? are poorly known 
quantities, separation of the theories does not appear possible based upon 
experiment. Numerical evaluation shows that the solution presented here 
lies between the Carstens and the Fukuta and Walter formulations, and 
even for a ~ A the solutions do not diverge by more than 7 percent. 
Thus, it appears that the differences in handling the physical assumptions 
which gave rise to the many different formulations in the literature are 
not important to the apphed problem. 
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APPENDIX 

EVALUATION OF THE INTEGRAL. A(a/A) 


If we write s for a/ A and t for r/a. then equation (8) in the text 
becomes 


„ /** e sin0 cosG 

l(S) = s # / K- 

0 0 (1 + 2t tx)se + t ) 


de dt 


and one integration yields 


A(S) = I e-*' ^ a. 


According to the Binomial Theorem, 




a . , 


t - i- + y 

' .2 L 2 ( 4 ) . 


. (3 - 2n) 1 

. (2 + 2n) ^2n + 1 


for Kt 


and, hence. 




s / V 3( 1) . 
3 I U 2(4) . 


. (3 - 2n) f ^-st *2n „-st * 

:~ (2 ^2n) I ® ^ I ^ ^ 
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■ /* ^ n4 Isl \ • 


or, with some combination of terms. 


A(8) -- I ^ '*•) 




n=0 

Now, if we let 


I lli ; : v irvinl (I «) 


(A-l) 


H8.n) =/ e-»' dt . / d, 

0 1 ^ 


(A-2) 


whenever n is non-negative and s is positive, (A-l) may be written as 


»y\_ s,/ l\ s ^ 3(1) . . . (3 - 2n) , 

A(s) = ' 5 <{< (^s, jj+ 5 2 5r?TTTTlTT*Iin 

n=0 


or, with the first two terms interchanged, 


A(8) = I ,j;(s.O) - I ,Ks,J) + I £ - 4 - B TT Th '~ (IV In ' ) 

n=l 


I ij^(s.O) - I ii)(s,l/2) + I < 1 .( 8 , 1) - ^ <Ks,2) + 


• • 


(A-3) 
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In view of (A-2), we see that 


0 < (s.n) < 


/q ^ ^in + T - Wn. ^ 5H n 


A e-st 

J 0 


.oo -St 

1 


dt +/~ e dt = r 


and that 


3i|>(s.n) 

3n 


I KAJ 

=/q t^ (2 In t) dt + /j 


-St 


(-2 In t) 


,2n + 1 


dt < 0 


Thus, i|;(s,n) decreases steadily to zero as n incrjases and, since the 
multipliers in (A-3) also decrease in magnitude, it is clear that (A-3) is 
a convergent alternating series. Moreover, each partial sum is an upper 
or lower bound for A(s) dep'^nding on whether the number of terms 
included is odd or even. For example: 

0 < A(s) < I i|»(s,0) 

1 'i^(s,0) - I i;;^s, < A(s) < I Ks,0) - I i|)^s, I ij;(s,l) 


and so on. 

The following well-known formulas may be verified by successive 
integration by parts for s ^ 0 and n > 0: 


/ t" dt 

0 


— 1 6 
s" * M 



OO 



-.-t 

e 


,n 


dt 




k 


s 


k! 
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and 


f . (-1)"'^ s"'^ L . . _ -a (-1)^ k!\ 

/j ^srrr ' (n - d! ® j|b s’' ^ / 


where 


00 00 -st 

EjU) = / 2j- at =/ V ■“ 

S 1 


We use these and equation (A-2) to obtain the alternative form 


<P(s,n) 


(2n)f 

„2n + 1 
s 




tijjn 8 

(2n) 


2n 


^E^(s) 


« 2nrl 

s 

k=0 


Lzll 




In particular, using A^(s) to denote the sum of the first n terms of 
equation (A- 3), we have 


Aj(s) = j \l>(a,0) * 


(A-4) 


A 2 (s) = Aj(s) - I ij;(s.l/2) = j - jg- 






. V s -s 

(8) - « C 


(A-5) 
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AgCs) = AgCs) + I ij^(8,l) 


. 1 
‘ 5 “ 


1 




ik - 


(A-6) 


(l*T *rs) *n) 


-s 


Substitution of a/X for s shows that equation (9) is the same as 
A 2 (s). By taking higher values for n, we can obtain approximations for 

A(s) to any desired degree of accuracy. For example, Ag(l) = 0.31933. 

Ag(l) = 0.31920, and we see that Ag(l) can be in error by no more than 

about 1/10 of 1 percent. Also, it can be shown that, if n is positive. 


Um 


n 


s-O ® = 


lim 

S-*-" 


n 


E^(s) = 0 


(A-7) 


and, consequently, the formvdas for A.(s) and A„(s) are sufficient to 
imply that ^ ^ 


lim 
s^ 0 


A(s) = 0 and A(s) 

s->“ 


1 

2 


(A-8) 


The numerical value of A^ts) is illustrated graphically in Figure 1. This 

third-order approximation to A(a/X) should be more than adequate for 
most practical applications of this droplet growth theory. 

As a by-product of the preceding derivation, we have obtained a 
Laplace transform of 


( 1 . - 1 


- t 


Generalizations are clearly possible for some exponents other than 3/2. 
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